Nuclear Rings in Galaxies — A Kinematic Perspective 
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ABSTRACT 

We combine DensePak integral field unit and TAURUS Fabry- Perot observations of 13 nuclear 
rings to show an interconnection between the kinematic properties of the rings and their resonant 
origin. The nuclear rings have regular and symmetric kinematics, and lack strong non-circular 
motions. This symmetry, coupled with a direct relationship between the position angles and 
ellipticities of the rings and those of their host galaxies, indicate the rings are in the same plane 
as the disc and are circular. From the rotation curves derived, we have estimated the compactness 
{v^ /r) up to the turnover radius, which is where the nuclear rings reside. We find that there is 
evidence of a correlation between compactness and ring width and size. Radially wide rings are 
less compact, and thus have lower mass concentration. The compactness increases as the ring 
width decreases. We also find that the nuclear ring size is dependent on the bar strength, with 
weaker bars allowing rings of any size to form. 
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1. Introduction 



Galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure - galaxies: 



Heller fc Shlosman I (|l994l Il996l): IPiner et al 



Nuclear rings in the central kiloparsecs of 
galaxies are usually observed as conglomera- 
tions of yo ung massive stars in d i stinct compact 
groupings (Buta fc Combes lll996tlKnapen Il2005t 



19951) : iKnaoen et al.l (119951): iButa fc Combes . 

1996l i: lRegan et all (|l997l ): [Benedict et al.l(|20o3 i) 

Some 20% of nuclear rings occur in non-ba rred 



Comeron et al.l 120101 ). The star forming rings 
are believed to form as a result of radial gas in- 
flow towards the central region, which stagnates 
near dynamical resonances. In barred galaxies 
the inflow originates in th e gravitational t orque s 
set up by the bar (e.g., lAthanassoula I (|l992| ): 



galaxies ( Knapen Il2005t IComeron et al.ll2010[ ) and 
in those cases the inflow and resonances may well 
be related to the presence of weak ovals, past in- 
teractions or mergers, or even strong spiral arms. 
A resonant origin of all nuclear rings is most plau- 
sible. 

The radial location of a nuclear ring can be 



described by the position of the inner Lind- 
blad resonancc(s) (ILRs) (ISimkin et ah' 'l980|: 
Shlosman ct al. 1989; Buta & Crocker 1993; Friedh 
19931: JHeher fc Shlosman I Il994 iKnapen et alT 



be an indicator of active star formation and cool 
_gas inflow from which the massive stars inherited 



19951: 

cally at 



Buta fc Combes 



19961 ). which is typi- 
the galaxy center al- 



1 kpc from 
though smaller a.nd larger nuclear rinRS exist 
(|Buta fc Crocker Ill993l : fComeron et al.ll2010[ ). In 
a non-axisymmetric potential, the nuclear ring 
forms near the ILR(s), where the inflowing molec- 
ular gas slows down, accumulates, and initiates 
star formation. For barred spiral galaxies, the 
ring can be formed inside the large scale galactic 
spiral structure in the vicinity of the inner Lind- 
blad resonances (ILRs), which occur at those radii 
where the bar pattern speed fib equals (fi — k/2), 
with Q the angular velocity and kappa the ra- 
dial epicyclic frequency. This approximation is 
not valid for strong bars where a non-linear anal- 
ysis must be performed, but in any case serves 
as a useful indication of how the location of the 
ring is related to the underlying galaxy dynamics 



fc Benm e low velocity dispersion. 

The existing kinematic studies provide insight 
into a small number of nuclear rings, but fail to 
reach more general conclusions applicable to the 
overall population. Apart from the differing char- 
acteristics of the data used in these studies, this is 
due mostly to the intrinsically small angular size 
of nuclear rings, and the dusty environments in 
which they exist, both of which make rings non- 
trivial to observe with most common instruments 
yielding kinematic data. In this paper we endeavor 
to start rectifying this situation by presenting a 
more homogeneous data set on thirteen nuclear 
rings. 

The nuclear rings studied spectroscopically in 
the current pa p er we re all part of the sample of 
Mazzuca et all (|2008[ ) (hereafter known as MOB), 



(JKnapen et al.lll995l ). 



Only a few in-depth observational studies of the 
kinematics of nuclear rings at high angular res- 
olution exist, but these have revealed some con- 
sistent kinemat i c char acteristics of nuclear rings. 



Benedict et al.l (|l996l ) showed that non-circular 
motions near the intersection of the ring and 
bar dust lanes correspond to radially inflowing 
molecular gas t hat is apparently feed i ng th e ring 
in NGC 4314. iRevnaud fc Downes I (Il998l) con- 
cluded for NGC 1530 the existence of weak ve- 
locity transitions at the ends of the bar which 
are inversely correlated with strong star forma- 
tion in th e same areas , includ ing the nuclear ring 



location. iRegan et al.l (|l997l ) further note an in- 
crease in the residual velocity alon g an axis per- 



pen dicular to the bar. [Knapen et al. (2000) (see 



also lKnapen et al.l (|l995[ )) and lZurita et al.l (|2004h 



used CO and Ha emission in the central region of 
NGC 4321 and NGC 1530, respectively, to reveal 
gas streaming i nward towa r ds the ring along the 
bar dust lanes. Ijogee et al.l (|2002l ) studied the fu- 
eling mechanisms of the ring in NGC 5248 to de- 
duce that large star forming clusters in the ring 
have been triggered by a bar-driven spiral density 



wave. 



More recently, lAllard et all (|2006l ) follow 



up on the ring in NGC 4321 to show relatively 
low H/? gas dispersion within the ring, which can 



who provided photometric insight into the strong 
star forming nature of nuclear rings. From their 
Ha observational study of 22 nuclear rings, they 
find that many of the rings exhibit a well-defined 
age distribution pattern. M08 conclude that this 
occurrence can be a result of the combination of 
bar-induced dynamics and gravitational instabil- 
ities which are occurring in the proximity of the 
rings. These results add credence to the idea that 
the mass infiow driven by the bar along the dust 
lanes to the ring contact points is a key require- 
ment for the fueling of active (i.e., star forming) 
nuclear rings (and might also yield a possible sig- 
nature in the form of non-circular motions) . Fur- 
ther evidence for this idea from detailed observa- 
tions of a handful of galaxi es has been presente d 



by, e.g.. lAllard et al.l (I2006D : iBoker et al.l (|2008f ). 



and lKnapen et al.l (J2010l) . 



We now use kinematic information within and 
just outside the nuclear rings in 13 galaxies to in- 
vestigate the possible relations between the optical 
morphology of the rings, the parameters derived 
from the rotation curve, and their resonant na- 
ture. We begin by presenting an overview of the 
observations and morphological properties of the 
Ha sample in Section 2. We discuss the basic data 
reduction method, which includes the construc- 
tion of the velocity field maps, in Section 3. In 
Section 4 we discuss the kinematic parameters as- 
sociated with the ring rotation curves, and present 
results pertaining to non-circular motions in and 



around the nuclear rings, including velocity resid- 
uals. Our detailed analysis of the results is dis- 
cussed in Section 5, with concluding remarks in 
Section 6. 

2. Sample and Observations 

The galaxies studied here (see Table [1]) form 
a subset of the sample presented in the imaging 
survey of M08 (which was based on B, I, and 
Ha imaging). The data in M08 were obtained to 
characterize the morphology of the nuclear rings, 
and to analyze the age distribution of the indi- 
vidual stellar clusters forming each ring using the 
Ha equivalent width. Results from this photo- 
metrically based survey reveal that the rings con- 
tain very young (1 Myr to 10 Myrs) hot massive 
stars, with some of the youngest clusters located in 
proximity to one or both of the intersection points 
of the bar dust lanes and the outer edge of the 
ring. The ring ellipticity and position angle ap- 
proximately match that of the respective galac- 
tic disk, indicating that the rings are in the same 
plane as the di s k and are circular (M08; see also 



Comeron et al 



(HqiS)). 

The kinematic sample was chosen from the 
imaging parent sample based on unambiguous de- 
tection of H II regions delineating the ring, a non- 
zero inclination (i), and a measurable ring posi- 
tion angle (0r). Rings too diffuse or small, or 
with no resolved Hii regions detected using the 
method described in MOB (i.e., weak signal and 
resolution), were rejected. With the Densepak in- 
strument we cannot detect any Hii regions for 
nuclear rings with angular sizes less than O'.'S in 
radius, however, we have previously unpublished 
high resolution Fabry-Perot Ha data of NGC 1300 
and NGC 6951 from the TAURUS instrument on 
the 4.2 m William Herschel Telescope, which we 
include in this sample. Similarly, we have used 
the TAURUS data for NGC 4321 from Knapen 
et al. (2000), a nuclear ring which meets our cri- 
teria but was added to the M08 sample after the 
DensePak observing runs. Lastly, although not in- 
cluded in the M08 sample, we observed two well 
known nuclear rings (NGC 2903 and NGC 3351) 
during the December 2003 run, which we also add 
to this sample. 

We obtained two-dimensional Ha velocity fields 
of the sample over two observing runs, in De- 



cember 2003 and April 2004, using the DensePak 
fiber-optic array and bench spectrograph on the 
3.5 m WIYINLj telescope. The DensePak integral 
field unit (IFU) spectrograph (JBarden fc Wade I 
19881 ) consists of 91 fibers bonded into a 7 x 13 



staggered rectangular grid that covers an area of 
30 X 45 arcsec of sky with center-to-center fiber 
(each fiber is 2.82 arcsec in diameter) spacings of 
3.75 arcsec. Four of the fibers are allocated as sky 
fibers, and are spaced around the grid roughly 1' 
from the center. The spectrograph was configured 
with the Bench camera, which uses a Tek 2048 
CCD (T2KC), with the 860 Zmm"! grating, in 
the second order, centered at 6575 A. We used a 
third order blocking filter (GG495). This arr ange- 
ment provided a final spectral coverage of 925 A, 
from 6085 A to 7010 A, at a resolution of 0.45 A per 
pixel. The spectral range included Ha as well as 
the [Nii]AA6548,6583and [Sii]AA6716,6731 lines. 
The effective instrumental resolution was different 
for the two runs, which is likely due to differences 
in the instrument setup. For the December 2003 
run, the instrumental resolution (FWHM) aver- 
aged O.97A (44 kms~^) around Ha, with the res- 
olution for the April 2004 run 1.31 A (60 kms~i). 

The array contains six dead fibers which de- 
creases the number of active fibers to 85. To avoid 
these areas, the projected image of each nuclear 
ring was placed in the upper half of the fiber ar- 
ray. Because of this method, many of the rings 
were close to the edge of the field of view (FOV) 
along the array minor axis, which limits our dis- 
cussion of the environment exterior to the rings. 
We used two to four pointings, each of 1800 s in 
duration, offset to fill the area of the dead fibers. 
This method also improved the sampling in some 
cases. For all exposures, the major axis of the 
array was aligned along a North-South direction. 

The Fab ry-Perot observations for NGC 4321 are 
detailed in iKnapen et al.l ( 20001) . with very sim- 
ilar procedures implemented for NGC 1300 and 
NGC 6951. The pixel size in the reduced Fabry- 
Perot cubes is 0.27 arcsec squared, and the sep- 
aration between velocity channels is 15.7 km s~^. 
Each cube has 55 planes, thus covering a total ve- 
locity range of just over 1000 km s^^. Continuum 



^ Joint facility of the University of Wisconsin-Madison, Indi- 
ana University, Yale University, and National Optical As- 
tronomy Observatories 



Table 1: Morphological and photometric parameters of DensePak sample 
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Run 
Date 
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(°) 


(°) 






(") 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



473 
1300 
1343 
1530 
2903 
3351 
4303 
4314 
4321 
5248 
5953 
6951 
7742 



SAB(r)0/a 
SB(rs)b 
SAB(s)b 
SB(rs)b 

SAB(rs)bc 
SB(r)b 

SAB(rs)bc 
SB(rs)a 

SAB(s)bc 

SAB(rs)bc 
SAa 

SAB(rs)bc 
SA(r)b 



-4.4 
ii/a 
-4.4 
-1.6 
-6.8 
-4.2 
-3.3 
-3.0 
n/a 
-3.3 
-6.0 
n/a 
+0.7 



+7.2 
n/a 
+7.8 
+11.3 
+8.3 
+6.1 
+7.5 
+6.3 
n/a 
+6.5 
+6.6 
n/a 
+7.5 



49 
40 
41 
48 
65 
34 
18 
26 
30 
46 
26 
45 
18 



154 
135 
60 
25 
8 

20 
88 
135 
170 
115 
172 
146 
133 



153 
106 
80 
8 
17 
13 



30 
110 
169 
170 



164 
102 
82 
122 
24 
112 
10 
135 
153 
137 

no bar 
85 

no bar 



0.37 
0.25 
0.25 
0.35 
0.4 
0.3 
0.14 
0.1 
0.12 
0.3 
0.1 
0.2 
0.05 



0.37 
0.34 
0.37 
0.29 
0.52 
0.32 
0.11 
0.11 
0.15 
0.28 
0.17 
0.17 




12.2 X 6.9 
4.1 X 3.1 
8.8 X 6.6 
6.8 X 4.9 
1.7 X 1.0 

2.7 X 1.9 
3.3 X 2.8 
6.6 X 5.9 

8.8 X 7.0 
6.6 X 4.6 
6.1 X 5.5 
4.6 x3.7 

9.9 X 9.4 



1.7 X 1.0 
0.3 xO.2 
1.2 X 0.9 
1.2 X 0.8 
0.06 X 0.04 
0.15 X 0.11 
0.2 X 0.2 
0.3 X 0.3 
0.7 xO.6 
0.7 X 0.5 
1.0 X 0.9 
0.5 xO.4 
1.0 X 1.0 



Dec03 
Sep98 
Dec03 
Dec03 
Dec03 
Dec03 
Apr04 
Apr04 
May95 
Apr04 
Apr04 
Sep93 
Dec03 



Notes: Morphological and photometric parameters for the observed sample. Galaxies are listed by NGC 
number in order of increasing RA (col. 1) with morphological type (col. 2) from de Vaucouleurs et al. (1991; 
hereafter RC3). Col. 3 lists the offset of the optical center of the Ha image from the center of the DensePak 
array; this is not apphcable to the TAURUS data for NGC 1300, NGC 4321, and NGC 6951. The nuclear 
ring inclination i (col. 4) and photometric ring position angle (j)^ (col. 5) have been derived from Ha imaging 
by M08). The disk position angle (pd (col. 6) is from the RC3 with a dash indicating that a ring is circular 
with no definable position angle. The bar position angle (p^, (col. 7) and photometric ring ellipticity Cr (col. 8) 
are from MOB. Disk ellipticity ed (col. 9) is from RC3. The radius of the ring, shown as the semi-major axis 
by the semi-minor axis (col. 10 and col. 11), was derived from the imaging data of M08 using the distance 
given there. T h e ring semi-major axis and ring ellipticity values are all consistent with those reported in 
IComeron et al.l ( 20101 ). taking into account that the methodology applied was slightly different. Run dates 
are in col. 12. 



subtraction was performed using the channels on 
either extreme end of this range. 

3. Data Reduction 

The raw CCD frames of the DensePak data 
were zero- and bias-corrected, flat-fielded, cosmic- 
ray cleaned, and combined using the standard 
IRAF NOAO packages (zerocombine, ccdproc, 
FLATCOMBINE, COSMICRAYS, and IMCOMBINE, re- 
spectively) . The combined images were wavelength- 
calibrated using the NOAO iraf data reduction 
package dohydra to extract the spectrum yielded 
by each fiber. A CuAr comparison source was 
used for wavelength calibration. The reduced 
two-dimensional spectral image contains 89 rows 
(comprising 85 data fibers and 4 sky fibers), where 
each row is one pixel in height and corresponds to 
a single fiber. 

We created GIPSMJ scripts to identify and 
then subtract the continuum and sky from the 
spectrum following t he method described in 



Bershadv et al.l (|2005| ). Due to pointing uncer- 
tainties inherent to the instrument, we established 
the measured pointing offsets (i.e., the shift from 
the center of the DensePak fiber array to the ring 
center) from the observations the mselves. W e 



follow ed the procedure outlined in ISwaters et al. 



( 20031 ) ■ in which the continuum levels measured 
from the spectra are compared with levels derived 
from a reference image with known coordinates. 
For our sample we used the /-band images from 
M08. By repeating this process for all positions 
near the nominal pointing, we created a grid of 
X^ values, in which the minimum y^ corresponds 
to the measured pointing offset. We estimate 
that the uncertainty on the derived pointing is 
approximately O.Sarcsec. Because the sky fibers 
contained no Ha emission from the target galaxy, 
we used the average of the four reserved DensePak 
sky fibers, each located 1 arcmin from the center 
of the array. 

To combine the pointings into a sparse velocity 
field, we first fitted Gaussians to the continuum- 
and sky-subtracted Ha emission line profiles, using 
a 3ct threshold to avoid spurious detections. The 
fitted velocities derived for each fiber were placed 
in a map at corresponding fiber positions. The 



observed radial velocities were corrected to helio- 
centric velocities. In the cases of multiple velocity 
measurements for the same pixel in the velocity 
field, we averaged the measured velocities. To aid 
in the visual representation of the sparse velocity 
field, we also constructed interpolated smoothed 
velocity field maps. We replaced each data point 
in the sparse velocity field with a Gaussian beam 
of 5.5 arcsec and weighted the overlapping sections 
between the points by the relative intensities of the 
overlapping Gaussians. The resulting smoothed 
velocity fields are shown in Figure [1] (left, panel 1). 
Because of interpolation effects these contiguous 
velocity fields may be uncertain, especially near 
the edges. Any quantitative interpretation comes 
from the sparse velocity maps, which are not in- 
terpolated. 

The Fabry-Perot data were wavelength and 
phase calibrated, with sky and continuum emis- 
sio n subtraction perform ed, as described in detail 
in iKnapen et al.l (2000|). The resulting reduced 



^ Groningen Image Processing 

|http://www.astro.rug.nl/gipsy| 



System, 



data cubes were used to produce Ha velocity field 
maps using GIPSY , and follo wing the pre scriptions 
given by Knap cn (|l997t ) and lKnapen et al. ( 200(j ). 
The resulting velocity fields are shown in Figure[5] 

4. Analysis 

4.1. Rotation Curves 

Rotation curves are powerful tools for revealing 
intrinsic kinematic characteristics of galaxies, in- 
cluding those of the nuclear rings they may host. 
In our sample, the range in radius spanned by 
the rotation curves is small when compared to the 
overall disk size, but this is a result of the small 
field of view of the instrument and of the sensitiv- 
ity limit. 

For each galaxy we derived the rotation curve 
(Figure [3|) by fitting tilted rings to the sparse 
velocity field th rough an iterative process (see 
iBegeman I (|l989l ) for procedural details) . We chose 
the incremental radius for each concentric ring to 
be 1.5 arcsec, which is approximately half of the 
width of a given fiber. We supplied the initial 
input values for the nuclear ring center, inclina- 
tion, and position angle, as measured from /-band 
images (M08). We derived the initial center posi- 
tion by taking the point where the /-band emission 
is the highest near the galaxy center. The input 
value for Vsys in the fit to the velocity field is from 






Fig. 1. — Results from the DensePak sample. From left to right, / band image from M08 with velocity field 
overlaid in contours, residual velocity field with contours, and Ha emission image from M08. The center of 
the galaxy is indicated with a cross in each panel. Contour levels in the left panel are spaced by 5 kms^^ 
with dotted contours representing velocities less than the systemic velocity and solid contours representing 
those velocities greater than the systemic velocity. Residual velocity contours and grayscales within the rings 
generally range from —10 to 10 kms~^ except for NGC 1530 whose contours range from —20 to 20kms"^, 
and NGC 473 and NGC 7742 which both have residual ranges from —5 to 5kms^^. N is up, E to the left. 
The size of the field shown is lOarcesc (RA) by 40arcsec (DEC) for all galaxies. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 2. — As Fig.l, now for the three galaxies observed with the TAURUS Fabry-Perot instrument. Panels 
are, from left to right, velocity field (grayscale and contours), result after subtracting a model velocity field 
based on the rotation curve from the observed velocity field, and Ha integrated intensity. The center of 
the galaxy is indicated with a cross. The field of view shown is 15, 35, and 15arcsec square for NGC 1300, 
NGC 4321, and NGC 6951, respectively. Contour levels for the first (leftmost) panel are for NGC 1300 from 
1472 to 1552kms"i (white contours) and from 1572 to 1672kms-i (black); for NGC 4321 from 1480 to 
1560 (white) and from 1580 to 1680 kms^^; and for NGC 6951 from 1300 to 1440 (white) and from 1460 
to 1600 km s"^ — all contours arc spaced by 20kms~^ and the grayscale range is the same as that spanned 
by the contours. In panel 2, contour levels for all galaxies are from —45 to Oknis^^ (white) and from 15 to 
45kms~^ in steps of 15kms~^; the grayscale range is from —50 to 50kms~^ for NGC 1300 and NGC 6951 
but from —30 to 30kms~^ for NGC 4321. For all three galaxies, the grayscale in panel 3 indicates the range 
of values in instrumental units. 



the NASA/IPAC extragalactic database (NED). 
One exception is NGC 5953, whose systemic veloc- 
ity was derived from the data due to a 26 kms~^ 
difference with the NED value. For the initial run, 
although we supplied input values, all parameters 
were left free. We first determined the center posi- 
tion and then fixed only that value for the second 
iteration to then determine the systemic velocity. 
Once the center and systemic velocity were fixed 
we determined the inclination and position angle 
(third and fourth iterations, respectively). With 
the four values then fixed, we fit the rotational 
velocity. We gave more weight to positions along 
the major axis, since those data points carry most 
of the information about rotational velocities. In 
most cases we used the inclination of the outer disk 
for the final fit. However, in the cases of low incli- 
nation (below 25°), we adopted the photometric 
orientation parameters from M08 as initial inputs 
to the routine since the fitting routine produced 
inclinations with high uncertainties. Even with 
such low inclinations the position angle was deter- 
mined accurately with a tilted ring fit. Because 
our data sets span a small range in radius, it was 
not feasible to make sophisticated fits to the ve- 
locity fields. 

The uncertainties on the velocities derived from 
the line fits range from 0.5kms^^ to 5kms^^. An 
additional source of uncertainty is due to the un- 
known location of the Halpha emission within the 
area seen by each fiber. The amplitude of this un- 
certainty depends on the velocity gradient across 
the face of the fiber. We estimated this gradient 
for all fibers using the derived rotation curves, and 
found that a value of 8 kms~^ is representative of 
a la uncertainty. 

Rotation curves for the Fabry-Perot data were 
similarly generated, but here, the smaller pixel size 
and full spatial sampling allowed us to fit the ro- 
tation curves with radial steps of 0.5arcsec from 
the center outward. 

4.2. Residual Velocity Fields 

We constructed residual velocity maps (Fig- 
ure m panel 2) for the velocity fields (both the 
sparse and interpolated versions for the DensePak 
data) by subtracting a smooth model velocity 
field, constructed from the rotation curve and the 
assumed orientation parameters. As with the in- 
terpolated velocity fields, the interpolated residual 



maps should only be used for general impressions. 
These maps, however, are useful because they give 
a better sense of patterns in the image. 

4-2.1. Within the rings 

We find that all of the nuclear rings in our sam- 
ple have low residual velocities (both approaching 
and receding) near the central ridge of the ring 
(ranging from ^^Skms"^ to '^20kms~^). The nu- 
clear rings in NGC 473 and NGC 7742 are not 
resolved enough to comment on here, as the asso- 
ciated residuals do not overcome the uncertainties 
associated with the measured velocity, wavelength 
calibration, and position of the Ha peak within a 
given fiber beam. Patches of approaching and re- 
ceding residual velocity peaks of '^lOkms^^ occur 
in the nuclear rings of NGC 1343, NGC 4314, and 
NGC 5953, and reach to near 20 kms~^ in the 
rings of NGC 1530 and NGC 5248. Since these 
values are close to the combined measurement un- 
certainty, they can provide no firm evidence for 
the existence and origins of non-circular motions 
(in NGC 4321 the situation is better, as discussed 
in detail by Knapen et al. 2000). 

A comparison of the kinematic and photometric 
centers of the sample reveals no strong asymme- 
tries. The degree to which the centers differ is pro- 
portional to the lopsidedness of the potential, and 
thus a g ood indicator of the amount of non-circular 
motion (JFranx et al.l 119941 ). After computing the 



difference between the kinematic and photometric 
center positions (see Table [2]), we see that the cen- 
ter positions agree to within the uncertainties of 
the photometric center measurements. This agree- 
ment, in combination with the small noncircular 
motions, suggests the nuclear rings are circular in 
nature. 

4.2.2. In the vicinity of the rings 

For some of the galaxies in our sample, our ob- 
servations cover part of the velocity field outside of 
the nuclear ring as well. This allows us to study 
the interaction between the bar and the nuclear 
ring in those cases. 

Sufficient, albeit limited. Ha radial coverage al- 
lows us to comment on the velocities of the barred 
galaxies near the outer edge of the nuclear ring, 
perpendicular and parallel to the bar major axis. 
In all cases, we can see the influence of the strong 
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Fig. 3. — Rotation curves with the radial range (i.e., width) of the nuclear rings shadowed. The turnover 
radius is denoted by the dashed vertical line. See Table [3] for ring widths. 
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Table 2: Kinematic parameters for the observed sample 



NGC 


14ys 


Offset 


'/'kin 


(/"kin — 0phot 


<i>& - '/'kin 


1 Ctrkin 


— Ctrphot 1 






X 


V 








X 


V 




(kms-1) 


(") 


(") 


(°) 


r) 


(°) 


(") 


(") 


(1) 


(2) 


( 


:3) 


(4) 


(5) 


(6) 




(7) 


473 


2135 


-5.0 


+7.5 


156 


2 


3 


0.6 


0.3 


1300 


505 


n/a 


n/a 


80 


55 


26 


n/a 


n/a 


1343 


2215 


-5.3 


+6.5 


7 


53 


73 


0.9 


1.3 


1530 


2450 


-1.6 


+11.0 


5 


20 


3 





0.3 


2903 


560 


-4.0 


+7.4 


5 


3 


12 


2.8 


0.9 


3351 


764 


-4.2 


+7.0 


30 


10 


17 





0.9 


4303 


1560 


-3.4 


+ 7.5 


117 


29 


— 


0.1 





4314 


979 


-4.0 


+5.5 


122 


13 


— 


1.0 


0.8 


4321 


1578 


n/a 


n/a 


150 


20 


120 


n/a 


n/a 


5248 


1149 


-2.9 


+6.3 


115 





5 


0.4 


0.2 


5953 


1991 


-6.0 


+6.7 


49 


57 


120 





0.1 


6951 


1423 


n/a 


n/a 


136 


10 


106 


n/a 


n/a 


7742 


1658 


+0.7 


+6.6 


131 


2 


- 





0.9 



Notes: NGC numbers of the target galaxies are listed in order of increasing RA (col. 1). Systemic velocity 
(col. 2), nuclear ring kinematic center position with respect to the center of the DensePak array (col. 3) 
and kinematic position angle (col. 4) were derived using the gipsy rotcur task. Typical uncertainties for 
l^ys are ±8kms~^, center position ±0.09arcsec, and position angle ±5°. The difference (absolute value) 
between the kinematic and photometric position angles (see Table[T]) is in col. 5, with the difference between 
kinematic and disk position angles in col. 6. The difference between kinematic and photometric centers is in 
col. 7. 
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bars as they create high residual velocities. How- 
ever, once they interact with the nuclear ring outer 
boundary, residual velocities sharply decrease and 
plateau within the nuclear rings. For NGC 4314, 
we see ordered velocity patterns increasing near 
the outer edge of the ring (both approaching and 
receding) close to both sides of the bar minor axis, 
where the bar dust lanes merge with the ring. 
Residuals approach 60kms~^ in both directions, 
with a strong velocity transition zone near the 
outer boundary of the nuclear ring, where values 
reach 90kms~^ (approaching) at the edge of the 
FOV. Although the spatial sampling is poor at the 
image edge, ou r observations do agre e with veloc- 
ity values from [Benedict et al.l (|l996l) . 

The residuals for NGC 5248 near the western 
side of the bar minor axis reach 40 kms~^, with a 
weaker-defined, although evident, transition zone 
as radii approach the outer ring boundary. These 
results are consistent with t he two-dim ensional 
CO velocity field presented bv ljogee et al. (2002 ). 
In both galaxies, good agreement with the litera- 
ture adds confidence that the results we see at the 
field edges are not artificial, even though our data 
are spatially sparse. 

The residual map of NGC 1530 provides us with 
the clearest picture of the interaction between the 
bar and the ring. In Fig. |4] we see large velocity 
excesses (~110kms~^) to the west of the nuclear 
ring, near its exterior edge. The high residuals are 
associated with the large scale bar. The strong 
residuals sharply decrease to 50km s~^ over a small 
transition zone near the outer boundary of the 
nuclear ring. Although we see no obvious resid- 
ual patterns near the bar major or minor axes of 
the nuclear ring, we note that the largest residual 
velocities (~120kms~^) across our sample occur 
in this galaxy, between the bar major and minor 
axes, on the west side just exterior to the ring 
outer boundary. A closer look in the vicinity radi- 
ally shows a very sharp transition at the location 
of the nuclear ring, where residual velocities drop 
to lOknis^^. Zurita et al. (2004) observe a nearly 
identical pattern in the corresponding Ha map of 
the circumnuclear environment. 

The TAURUS data of NGC 4321 used here were 
analyz ed and discussed in detail bv lKnapen et al 
( 20001 ) ■ In the residual velocity map, they high- 
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Fig. 4.— Velocity curve for NGC 1530 at radii 
within (top plot) and exterior to (bottom plot) 
the nuclear ring. The radial coverage within the 
nuclear ring is further divided to show velocities 
along the inner and outer halves of the ring. Out- 
side of the nuclear ring, velocities from the nu- 
cleus to the ring inner edge are plotted, as well 
as those from the ring outer edge to the radial 
limits of the FOV. The velocity is plotted with re- 
spect to the systemic velocity, normalized by the 
rotation velocity, and adjusted by removing the 
inclination in order to allow comparison with the 
curve representing circular motion. The bar ma- 
jor axis (solid lines) and minor axis (dotted lines) 
are drawn for context. Note that there is a phase 
shift of 90°, plus the difference between the kine- 
matic position angle of the velocity data points 
and the photometric position angle used to cre- 
ate the model curve, to allow the cross-over point 
to occur at 180°. The area azimuthally located b 
etween the bar major and minor axes shows the 
largest velocities of the sample, with deviations 
reaching 110 km s^^. 



lighted the contributions from the inner part of 
the bar and the spiral arm fragment coming into 
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the nuclear (pseudo-)ring from the main part of 
the bar. These result in a broadly symmetric pat- 
tern of positive and negative residual velocities, 
mapping out the non-circular motions due to bar 
and spiral streaming motions. 

5. Results 

5.1. Circular Nuclear Rings 

For most of the sample, we find that the nuclear 
rings are nearly circular and in the same plane as 
the disk. The shape and orientation with respect 
to the host disk were deduced by comparing the 
position angle (of the nuclear ring, 0^ and disk, 
(pd) combined with the ellipticity of each nuclear 
ring, er to that of its host disk, Cd- We used the 
photometric position angles and ellipticities of the 
nuclear rings from M08. The kinematic position 
angles of the nuclear rings were extracted from 
the rotation curve generation. In Figure [SJ we 
compare the various position angles and see that 
in most cases the photometric and kinematic po- 
sition angles are similar. The same can be said 
for both the photometric and kinematic position 
angles with respect to the galactic disk position 
angle. This finding, coupled with the result from 
M08 of a similar relationship with respect to the 
disk and nuclear ring ellipticities (fourth panel), 
corroborate s our idea that nuclear rings are cir- 
cular and oriented as their host disk. 

We note three cases that are significantly above 
the scatter and thus do not meet our criteria 
of a round in-plane nuclear ring: NGC 1300, 
NGC 1343, and NGC 5953. The ring in NGC 1300 
is not fully sampled azimuthally (see Figure [2|) , 
which adds a high level of uncertainty for the pho- 
tometric position angle and ellipticity. NGC 5953 
is nearly face on and morphologically patchy and 
round along the line of sight, a combination that 
creates high position angle uncertainty. 

The residual maps (Fig. [1] and Fig. [U panel 2) 
also indicate that within the uncertainties of 
our observations (see Section 4.1), the nuclear 
rings are circular. There may be noncircu- 
larites at a lower level which may be related 
to the local underlying po t ential (as noted for 
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Fig. 5. — Comparisons for ring kinematic vs. ring 
photometric (top), ring photometric vs. disk pho- 
tometric (middle), and ring kinematic vs. disk 
photometric (bottom) position angles. Ellipticity 
comparison of the ring and disk is also plotted. 
The disk position angle is not given in the RC3 for 
NGC 4303, NGC 4314 and NGC 7742. The mean 
error on the disk position angle is unknown, but 
the error increases as the ratio of isophotal major 
to minor axes decreases (logi?25)- We assigned a 
5° error for the smallest ratio in the sample and 
proportionally computed the other errors. Errors 
are smaller than the size of the graphic icons of 
each galaxy. 
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Table 3: Nuclear ring and bar properties. 



NGC 


Nuclear 


RelRing 


Ring 


Q. 


Compactness 


Rt 


Disk 




Activity 


Size 


Width 










sma 








(kpc) 




(10'*(kms-i)2kpc-i) 


(kpc) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


473 


- 


0.20 


0.62 


0.12 




1.1 


1.4 


7.4 


1300 


- 


0.02 


0.13 


0.54 




16.5 


0.31 


18.2 


1343 


— 


0.11 


0.31 


0.15 




2.3 


0.3 


11.1 


1530 


— 


0.05 


0.62 


0.73 




8.2 


1.4 


24.5 


2903 


Hii 


0.08 


0.01 


0.27 




11.7 


0.1 


16.3 


3351 


Hii 


0.04 


0.09 


0.22 




36.1 


0.2 


4.9 


4303 


Hii 


0.01 


0.20 


0.28 




23.8 


.2 


14.4 


4314 


LINER 2 


0.05 


0.11 


0.43 




8.6 


0.3 


5.9 


4321 


Transition 2 


0.04 


0.38 


0.22 




15.8 


0.2 


21.7 


5248 


Hii 


0.03 


0.49 


0.10 




4.6 


0.5 


20.5 


5953 


LINER 2, Sy2 


0.13 


0.90 


0.1 




9.8 


0.35 


7.7 


6951 


Sy 2 


0.03 


0.23 


0.28 




9.9 


0.25 


13.8 


7742 


Transition 2 


0.18 


0.63 


0.06 




5.8 


0.32 


5.5 



Notes: Identification (col. 1); presence and type of nuclear activity (col. 2, from Ho et al. 1997b when available, 
NGC 1530 and NGC 5953 are from NED); relative nuclear ring size (col. 3) defined as ring semi-major axis (see 
Table[T] col 11) divided by disk semi-major axis (half of the host major disk diameter), Do from NED (col 9); ring 
width (col. 4) defined as the difference between the isophotal inner ring boundary and outer ring boundary, divided 
by the disk semi-major axis; non-axisymmetric torque parameter (col. 5), from Comeron et al. 2009; when no Qb 
value exists, we adopt a known Qg value instead: NGC 1530 (jBlock et al.l 120041 ) and NGC 5953 (H. Salo, private 
communication); rotation curve compactness (col. 6); rotation curve turnover radius (col. 7); We estimate that on 
average the uncertainty in the compactness is ~75kms~^kpc. Uncertainties in the various size measurements involved 
were also added in quadrature to yield typical uncertainties in the relative ring size and width of ±0.01. Uncertainties 
associated with the non-axisymmetric torque parameter, Qg, were adopted from Comeron et al. (2009) to be ±0.1. 
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5.2. Mass Concentration 
Ring Size 



and Nuclear 



An indication of the central mass concentra- 
tion can be o b tained from the compactness, e.g., 
de Blok et al.l (|l996l ). We define compactness as 
V^/Rt, where V is the rotational velocity differ- 
ence between the radial location of the turnover 
radius, Rt, and the origin (0,0). The turnover ra- 
dius is defined as the point in the rotation curve 
halfway between the initial steeply rising compo- 
nent and the flat(ter) segment. For the few galax- 
ies with very steeply rising rotation curves, such 
as NGC 2903 and NGC 3351, the turnover ra- 
dius cannot be measured accurately, and hence the 
compactness is uncertain. In these cases, the com- 
pactness should be considered a lower limit. 

Table [3] lists the compactness and turnover lo- 
cations. As expected, the compactness decreases 
with increasing turnover radius (not shown graph- 
ically, but the values are in the Table. In linear 
theory, which is not strictly applicable here but 
yields useful approximations, the turnover radius 
in the rotation curve sets the location of the ILRs, 
and they in turn determine where inflowing gas 
piles up, which ultimately forms the nuclear ring 



1996h . In- 



(jKnapen et al.lll995t iButa fc Combes 
deed, the location of the nuclear rings correlates 
well with the turnover radius (see Figure[3]and the 
left panel of Figure [6]), which we can identify with 



the phenomenon identified above. 

In panel 1 of Fig. [7] we compare the relative 
ring size, that is the nuclear ring semi-major axis 
divided by the semi-major axis of the host galaxy, 
to the compactness parameter as derived from the 
rotation curve. We find that the relative ring size 
increases as the compactness decreases. 

The relative ring size is compared to the non- 
axisymmetric galactic torque, Qg, in panel four of 
Fi g. E This torque pa r amete r, originally defined 
by iCombes fc Sanders I (|l98l[ ) to characterize bar 
strength, quantifies the overall influence of galactic 
non-axisymmetries. Specifically, Qg is the maxi- 
mum tangential force at a given radius divided by 
the radial force. For all galaxies except two, we 
adopt the Qg values from Comeron et al. (2010), 
as noted in Table[3l since they applied a consistent 
computation method for all of the nuclear rings 
seen in our sample. We confirm (Fig. [7]) that the 
highest torque values occur in those galaxies with 
the smaller rings relative to the size of their host 
disk (Knapen 2005; Comeron et al. 2010). As the 
torque weakens below Qg = 0.4 any size ring can 
form, but higher torque values ("stronger bars") 
progressively limit the nuclear ring size. 

The compactness and the non-axisymmetric 
torque Qg do not correlate (panel 3 of Fig. [T]). 
Also, the turnover and Qg do not correlate (not 
shown graphically) . This may be because the non- 
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axisymmetric torque is a parameter driven more 
by the azimuthal than the radial distribution of 
mass, whereas the opposite is the case for the 
turnover radius and compactness. 

5.3. Nuclear Ring Width Comparisons 

The width of a nuclear ring is an important pa- 
rameter because it very much defines how well a 
ring is visible in an Ha image, and how it is seen 
in contrast with the background emission from the 
host galaxy. Using the ring parameters we have 
now determined, it is interesting to see whether 
the ring width is driven by the physical properties 
of the host galaxy. So, as a first step, we define the 
ring width as the difference between the inner and 
outer ring boundaries, as defined from the emis- 
sion in Ha, and then plot the width of the nuclear 
ring against the relative ring size (right panel of 
Fig. in]). We see a trend that larger nuclear rings 
are wider, which might not have been expected a 
priori. 

We then compare the compactness to the nu- 
clear ring radial width (panel 2 of Fig. [7]). There 
is some evidence for a trend of compactness with 
ring width, where a nuclear ring is wider as the 
compactness decreases. We can thus conclude that 
increased compactness leads to both smaller and 
narrower nuclear rings, which is in accord with the 
expectations from linear theory and the location 
of a pair of ILRs near the turnover radius of a 
rotation curve. This clear link between the kine- 
matics of the galaxy and the morphological pa- 
rameters of the nuclear rings is thus further direct 
evidence for a resonant origin of the nuclear rings, 
deeply linked to the underlying structure of the 
host galaxy. This evidence based on kinematics 
supplements earlier evidence based almost exclu- 
sively on morphology, as reviewed and extended 
by Comeron et al. (2010). 

6. Concluding Remarks 

We combined DensePak integral field unit and 
TAURUS Fabry-Perot observations for a sample of 
13 galaxies that contain nuclear rings. Based on 
this sample, we note the following new findings: 

• Nuclear rings are intrinsically circular and 
unaffected by the local non- axisymmetric en- 
vironment that can exist from active star for- 



mation within the rings; they are infiuenced 
by global phenomenon that can affect the 
potential, such as a bar. 

• The compactness, derived from the rotation 
curves, decreases as the ring width increases. 

• The compactness also decreases as the rela- 
tive ring size increases 

• As the relative nuclear ring size increases so 
does the width of the ring 

• Strong bars with a Qb of 0.4 or greater can 
only be accompanied by small and thin nu- 
clear rings. For smaller torques, the ring size 
can vary. 

We have also reaffirmed several standing find- 
ings: 

• True to resonance theory, the location of 
the nuclear rings correlates well with the 
turnover radius. 

• The highest torque values occur in those 
galaxies with the smaller rings relative to the 
size of their host disk. 

In those cases with sufficient radial coverage, 
we agree that 

• For NGC 4314 residual velocities sharply de- 
crease and plateau within the nuclear rings. 
Ordered velocity patterns increase near the 
outer edge of the ring (both approaching and 
receding) close to both sides of the bar minor 



• The residual velocities for NGC 5248 near 
the western side of the bar minor axis reach 
40kms~^, with a weaker-defined, although 
evident, transition zone as radii approach 
the outer ring boundary. 

• In the case of NGC 1530 large velocity ex- 
cesses exist to the west of the nuclear ring, 
near its exterior edge and then sharply de- 
crease over a small transition zone near the 
outer boundary of the nuclear ring. 



The Fabry-Perot observations were made with 
the William Herschel Telescope operated on the 
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island of La Palma by the Newton Group of Tele- 
scopes in the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de 
Canarias. 

The Integral Field Unit observations were made 
with the WIYN Observatory which is owned and 
operated by the WIYN Consortium, consisting of 
the University of Wisconsin, Indiana University, 
Yale University, and the National Optical Astron- 
omy Observatory (NOAO). The WIYN is part of 
the Kitt Peak National Observatory located in 
Tucson Arizona. 
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